Scaling of hysteresis loops at phase transitions into a quasiabsorbing state 
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Models undergoing a phase transition to an absorbing state weakly broken by the addition of 
a very low spontaneous nucleation rate are shown to exhibit hysteresis loops whose width AA 
depends algebraically on the ramp rate r. Analytical arguments and numerical simulations show 
that AA ~ r K with k = l/(/3' + l), where /3' is the critical exponent governing the survival probability 
of a seed near threshold. These results explain similar hysteresis scaling observed before in liquid 
crystal convection experiments. This phenomenon is conjectured to occur in a variety of other 
experimental systems. 
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Directed percolation (DP) is an archetypical model of 
phase transitions into an absorbing state, i.e. a state from 
which a system can never escape. A vast literature of 
theoretical and numerical studies has enlarged the range 
of phenomena in the DP universality class [l[, refining 
conditions for this prominent critical behavior, known 
as DP conjecture [if, d, Q. Experimentally, the author 
and coworkers recently found that electrohydrodynamic 
convection of nematic liquid crystal shows the scaling 
behavior of DP at the transition between two turbulent 
states (DSM1-DSM2) 0|. Applying voltages V closely 
above the threshold, spatiotcmporal intcrmittency (STI) 
occurs, in which DSM2 patches move around in a DSM1 
background. As conjectured early by Pomeau [f|, this 
STI was unambiguously mapped onto DP with DSM1 
playing the role of the absorbing state. This constituted 
a clear experimental realization of a DP-class absorbing 
phase transition. 

On the other hand, Kai et al. reported in 1989 hys- 
teresis phenomena around this DSM1-DSM2 transition 
0| . Measuring the global light transmittance through 
the sample, increasing or decreasing the applied voltage 
V at a rate r, they found hysteresis loops of width AV 
scaling roughly like AV ~ r K with k ~ 0.5-0.6 [1,0|- In 
particular, these loops disappear in the small-r limit, and 
it has been discussed whether the transition corresponds 
to a supercritical bifurcation or a subcritical one. This is 
in apparent contradiction with DSM1 being an absorb- 
ing state, since then one expects infinitely wide hysteresis 
loops. It is shown here that the scaling of hysteresis loops 
is in fact in full agreement with the DP framework in 
which the DSM1 state is only quasi-absorbing, i.e. with 
the existence of a small residual probability for sponta- 
neous nucleation of DSM2 patches either in the bulk or 
at the boundaries. 

As a first illustration, a probabilistic cellular automa- 
ton (PCA) version of the contact process (CP) [HQ is in- 



troduced, in which an extra, small probability h to create 
an active site spontaneously anywhere is added. Consider 
a two-dimensional (2D) square lattice of size L x L and 
assign a variable s^j to each lattice point, encoding its 
local state, either inactive (absorbing, Sjj = 0) or active 
{si.j = 1). Indices i and j denote Cartesian coordinates. 
The time evolution is as follows: randomly choose one 
site and stochastically flip it with probabilities 



Pij(0 -► 1) = -^-(sj-ij -I- s i+ ij 

Pi j(l -^0)=P2, 



(1) 
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where p x = A/(A + 1) and p 2 = 1/(A + 1). The two 
terms in the first equation account for contamination by 
neighbors and spontaneous nucleation of active sites, re- 
spectively. Periodic boundary conditions Sij — Si+L,j — 
S i.j+L are used throughout, and a time step (or Monte 
Carlo step, MCS) consists of L 2 flipping attempts. The 
h = case is known as the PCA version of the orig- 
inal (2+l)D CP, which shows a DP-class transition at 
A c = 1.64877(3) [9| (the number in parentheses denotes 
the uncertainty in the last figure). In the present study, 
L = 256 and hi = hL 2 = 10~ 2 . Although, strictly speak- 
ing, even rare nucleation events wipe out the absorbing 
phase transition, in practice a significantly low nucleation 
rate allows us to observe the underlying critical behav- 
ior as we shall see in this study. The nucleation rate h 
theoretically corresponds to an external field [lOj . so a 
weak-field case is dealt with here. 

The model behaves similarly to the turbulence of liq- 
uid crystals in many aspects. For instance A > A c and 
initial conditions of Si.j — everywhere lead to a nu- 
cleus growth after sufficient time has passed, which faith- 
fully reproduces experiments. In particular, the model 
exhibits hysteresis as shown in Fig. [lja) an d Movie SI 
[ill ] when A is increased from A < A c to A > A c at a 
constant ramp rate r and then decreased at the same 
speed. The hysteretic process can be decomposed into 
three stages as indicated in the bottom of Fig. [2] Let 
us start from the uniformly inactive state and increase 
A. First, active clusters do not emerge even for A > A c 



2 




FIG. 1: Hysteresis observed in simulations and experiments, (a) Hysteresis of (2+l)D CP with h! = 10 -2 . Black regions 
denote active sites. The control parameter A is increased and then decreased in the range of 1 < A < 3.4 at the ramp rate 
of r — 0.001 MCS -1 . The critical point for the model without nucleation is A c = 1.64877(3) 9|]. The arrows after the values 
of A denote whether they are increasing or decreasing, (b) Hysteresis in the electrohydrodynamic convection, where the same 
cell as in Ref. is used. The control parameter, applied voltage V, is ramped in the range of 22 V < V < 75 V at the rate 
of r = 1.71 V/s with fixed frequency of 250 Hz. The critical voltage is V c w 35 V [J. Darker regions correspond to DSM2, the 
active state. Note that the global intensity and contrast are adjusted for the sake of clarity, and that DSM1 and DSM2 coexist 
in the two images at the lower left. 
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FIG. 2: (Color online) Typical hysteresis loops for four differ- 
ent ramp rates r in (2+l)D CP with ti = 10~ 2 . Note that the 
ratios of the four values of r are chosen to be approximately 
the same as in Fig. 1 of Ref. to allow the comparison (see 
also Note [HI)- The hysteretic process can be decomposed 
into three stages as indicated in the bottom figure. 



due to the very low nucleation rate (1st stage). How- 
ever, once a spontaneous nucleation occurs, the active 
nucleus grows and finally covers the whole system be- 
cause of A > A c (2nd stage). The density of active sites, 



p, saturates at the steady state value p 

steady 

(A). On the 

other hand, when A is decreased, the number of active 
sites decreases gradually and homogeneously contrary to 
the growing process, approximately following Psteady(A) 
(3rd stage). This strikingly resembles what is observed in 
the liquid crystal experiments [Fig. QJb) , Movie S2 [ll[ , 
Refs. [a 03]. Note that the observed hysteresis both in 
the experiments and in the simulations is not a station- 
ary property of the system, as would imply a first order 
transition, but rather a dynamical effect owing to the 
sweep of the parameter. 

The dependence on the ramp rate r is shown in Fig. 
[51 which is again very similar to the corresponding ex- 
periments @, [l2| • The widths of the hysteresis loops AA 
and AA*, defined as in Fig. [2] clearly exhibit the power 
law dependence AA, AA* ~ r K [Fig. [3] (disks and trian- 
gles)], with k = 0.61(1) for AA and « = 0.56(3) for AA*. 
Here the ranges of error correspond to 95% confidence 
intervals in the sense of Student's t. They are in good 
agreement with the experimental value re = 0.5-0.6 @, 0]- 

Besides the agreement between the simulations and the 
experiments, the exponent re can also be derived only by 
assuming DP criticality with a very low probability for 
spontaneous nucleation. For absorbing phase transitions, 
the probability P^ with which an active site survives 
forever grows algebraically as Poo ~ e 13 for e = A — 
A c > 0, where f3' constitutes one of the critical exponents 
characterizing these transitions. (Note that for the DP 
class the so-called "rapidity" symmetry implies /?' = 
1, 13], where j3 is the critical exponent corresponding to 
the stationary active site density p s toady) Suppose e is 
increased linearly as e(t) = rt and a nucleus appears and 
grows at time t — T, and assume that the ramp rate r is 
so slow that the finite-time survival probability converges 
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, k = 0,61(1) 

k = 0.56(3) 
k = 0.64(5) 
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FIG. 3: (Color online) Widths of the loops AA (disk), AA* 
(triangle), and A A* (square) with respect to the ramp rate 
r, in the case of (2+l)D CP with h' = 1(T 2 . The sym- 
bols and errorbars indicate means and standard deviations, 
respectively, of 50 independent runs. Dashed curves denote 
the results of the fitting to the power law AA, AA*, AA* ~ r K . 
The inset shows the same data in logarithmic scales. 



to Poo before the control parameter significantly changes, 
the following relation then approximately holds: 



1 



tiP oc (e(t))dt~tir' 3 'TP'+\ 



and thus the width of the hysteresis is 



AA* = rT 



(2) 



(3) 



It gives the exponent for the hysteresis as k = l/(/3'+l) = 
0.632(2) for the (2+l)D DP [14]. Of course the assumed 
nucleation process is stochastic, so that, strictly, one 
should deal with the average width (AA*) based on the 
probabilistic distribution. This more rigorous approach 
is also straightforward. With Po(t) being the probabil- 
ity that a nucleus does not appear and grow until time 
t, the probability that such a nucleation first occurs be- 
tween time t and t + At is written as 



-dP {t)=P {t)-h'P oo {e{t))dt 

( Ch'r ' 



exp 



t fr+1 dt, (4) 



where C is defined by ?«, = Ce@ . This gives the average 
of the hysteresis width as 



(AA t *) = r 



0> + 1 



(/3' + l)rl 1/(/3 ' +1) 



Ch' 



(5) 



which confirms Eq. ([3]). Note that the standard deviation 
also obeys the same power law (with a different coeffi- 
cient), since the stochastic process at play is essentially 
Poissonian. 



TABLE I: Hysteresis exponent re for several models. 



Model" 



Exponent re for 



~AA^l/(/3' + l) 



(2+I)D CP (PCA) 0.61(1) 0.56(3) 0.64(5) 0.632(2)^" 

(2+l)D CP 0.61(2) 0.61(4) 0.65(7) 0.632(2)'' 

(1+1)D CP (PCA) 0.69(1) 0.73(3) 0.81(4) 0.783 6 

(1+1)D site DP C 0.68(2) 0.71(5) 0.82(7) 0.783 6 

(2+l)D voter-like de 0.465(14) 0.460(17) 0.47(4) 0.5 



"System sizes and nucleation rates are set to L = 4096 and h' = 
1CT 4 for (1+1)D, and to L = 256 and h 1 = 10~ 2 for (2+l)D. 

'Values of the DP exponent /3' are from Ref. [l5ll for (1+1)D and 
from Refs. [JJ] for (2+l)D. 

Simulations are performed in the Domany-Kinzel lattice. 

d Kinetic Ising model with spin-flip probability p s n is considered, 
where s = ±1 and H € {—4, —2, 0, 2, 4} denote a spin and its local 
field, respectively. p_2 is swept here with the other parameters 



fixed at p4 



10- 2 /L z ,p 2 = 0.17, po = 0.5,p_ 4 = 0.68. This 



model shows a transition in the voter universality class |16|| . 

E Hysteresis is measured in terms of the density of interfaces (i.e., 
the fraction of H — pairs) instead of the active site density (i.e., 
magnetization), since the former characterizes the voter class better 
[lal and shows faster relaxation. 



The derived value of k = 0.632(2) is slightly larger 
than the numerical and experimental values. This stems 
from the use of different definitions for the lower bound 
of the loop: the merging point of the two curves defining 
the loop is used for the experiments and simulations (AA 
and AA*), whereas, theoretically, the exact critical point 
A c is used to define the lower bound. Adopting the latter 
definition for the simulations (AA* in Fig. [2]), k — 0.64(5) 
is obtained (Fig. [3]), which is now in close agreement with 
the theoretical value. Thus the picture based on DP with 
a weakly broken absorbing state quantitatively explains 
the observed hysteresis. The above derivation also indi- 
cates that the scaling of hysteresis loops is seen for such 
values of h! and r that nucleations, including those with a 
short lifetime, occur several times within the range where 
the scaling P^ ~ e 13 holds. 

Given that the only assumption was criticality of an 
absorbing transition together with very rare spontaneous 
nucleations, the observed scaling of hysteresis with the 
exponent k = l/(/3' + 1) is expected to be found uni- 
versally in systems that exhibit quasi-absorbing transi- 
tions. This is confirmed by performing simulations in 
different dimensions and for different models and univer- 
sality classes [Table Q] and Fig. Q] (a)(b)]. In all cases, the 
measured k values for AAjf are in good agreement with 
those derived by Eq. ^ . The loop scaling is also robust 
to situations when the ramp rate r and/or the nucleation 
probability h! vary with time or the control parameter. 
As long as they are nonzero and analytic at criticality, 
this gives only higher order corrections to Eq. © and 
does not affect the final result when r — > 0. Even if this 
condition is not satisfied, the corrected form of Eq. © 
can be calculated, for example in the case of nonlinear 
ramping e(t) = r't a ; the hysteresis exponent becomes 
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FIG. 4: (Color online) Hysteresis scaling for (1+1)D CP (a), 
(2+l)D voter-like(b), (2+l)D CP with quadratic ramping 
e(t) = r't 2 (c), and (2+l)D CP with square root ramping 
e(t) = r't 1 / 2 (d). The estimates of the exponent k for AAj 
agree with the theoretical values, namely, 0.783 (a), 0.5 (b), 
0.462(2) (c), and 0.774(2) (d). 



then k = X/(a/3' + 1), which is numerically confirmed 
[Fig.H(c)(d)]. 

Some experimental systems expected to belong to the 
DP class seem to lack strictly absorbing states due to 
residual nucleations [l], G3 • This suggests that the same 
hysteresis may be observed in such systems, for example 



with different alignments or at other transitions in the 
electrohydrodynamic convection (T5 , . A much more 
intriguing candidate can be found in the field of quan- 
tum turbulence (20| . Recently, a number of experimen- 
tal studies on transitions to turbulence in superfluid He 
have reported hysteresis [HI, [H, US H3] and temporal in- 
termittency in local state of turbulence [13, US HH HH . 
The existence of a (quasi-)absorbing state is also ex- 
pected due to the quantum topological constraint. All of 
these facts suggest that an absorbing transition to STI 
may take place in this superfluid system. Although it 
seems technically difficult to examine conventional crit- 
ical phenomena of absorbing transitions directly there, 
scaling of hysteresis loops may be more easily accessi- 
ble and would allow to decide about the corresponding 
universality class. 

In conclusion, the hysteresis loop scaling experimen- 
tally observed before at the DSM1-DSM2 transition of 
liquid cristal convection was explained by assuming DP 
dynamics with very rare spontaneous nucleations. This 
implies that DSM1 is probably only quasi-absorbing in 
the liquid crystal system. Moreover, scaling of hysteresis 
loops AA ~ r K with k = + 1) was demonstrated to 

be able to decide the universality class of transitions into 
a quasi-absorbing state. These results may also be used 
to analyze critical phenomena in systems where measur- 
able quantities are so limited that usual approaches to 
absorbing phase transitions cannot be adopted, such as 
in superfluid turbulence. 
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